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Abstract: We construct an entangled photon polarimeter capable of 
monitoring a two-qubit quantum state in real time. Using this polarimeter, 
we record a nine frames-per- second video of a two-photon state's transition 
from separability to entanglement. 
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1. Introduction 

Photonic entanglement is a fundamental resource for quantum information processing and 
quantum communications Q • Engineering suitable entanglement sources for a particular ap- 
plication, or integrating those sources into a larger system, however, can be a challenging 
experimental task. Generating high-quality entanglement requires protecting against or com- 
pensating for decoherence, single-qubit rotations, and partial projections. For both free-space 



E El SI 3 and fiber/waveguide-based entanglement sources (6llU[8][9l[T0]], this means com- 
pensating for any polarization rotations or decohering effects which may occur in transit to a 
destined application. In addition to the aforementioned static effects, it is necessary to test the 
source's stability in the face of real-time system perturbations such as atmospheric turbulence 
or fiber breathing owing to environmental fluctuations. At present, the best available technique 
for measuring two-qubit entangled states is quantum state tomography (HI El, a procedure 
which can provide a precise reconstruction of the quantum state, but which generally requires 
5-30 minutes to complete. This long measurement time can make debugging systematic exper- 
imental problems — particularly those with short timescales — challenging, if not impossible. 

The field of classical optical communications has faced similar problems when transmitting 
polarized light over long distances. A polarimeter is a common tool which is used to debug 
unwanted polarization rotations or depolarization effects (i.e., polarization decoherence). A 
polarimeter actively monitors the polarization state of a classical optical field, providing an 
experimenter with a real-time picture of the optical field's Stokes vector (i.e., its polarization 
state). Similarly, an entangled photon polarimeter — a measurement device capable of perform- 
ing quantum tomographies and displaying the reconstructed two-qubit states in real time — 
would be a valuable tool for optimizing and deploying entangled photon sources. 

In this paper we present the first experimental implementation of an entangled photon po- 
larimeter, which is capable of displaying nine reconstructed density matrices per second via 
complete quantum state tomographies. This represents a speed improvement of 2-3 orders of 
magnitude over the best quantum state tomography systems currently in use in laboratories 
around the world. Using this new tool, we record the first live video — at 9 frames-per- second 
(fps) — of a two-photon quantum state's transition from separability to entanglement. 

2. Two-Qubit Polarimetry 

Two-qubit polarimetry is a specific example of two-qubit quantum state tomography, a proce- 
dure for reconstructing an unknown quantum state from a series of measurements (generally 
either 9 or 36 coincidence measurements performed using two single-photon detectors per qubit 
lfT2l ). each performed on an ensemble of identical copies of the unknown state. Three key pa- 
rameters can be used to characterize any experimental apparatus for quantum state tomography: 
the time required to complete the state reconstruction procedure and the accuracy and precision 
with which the reconstructed density matrix represents the unknown quantum state. 

The time required to complete a quantum state tomography, T, is dependent on the number of 
two-qubit measurement settings taken per reconstruction, M; the time per measurement setting, 
T m ; the time to switch between measurement settings, T s \ and the time necessary to numerically 
reconstruct the unknown density matrix from an analysis of the measurement results, x a \ 

r=Mx(T m + T,) + V (1) 

The accuracy and precision of a tomography are closely related, both indicating how closely 
the reconstructed density matrix, p, matches the "true" unknown density matrix, Pideal- The 
"accuracy" of a tomographic reconstruction measures error due to systematic effects, such as 
improperly performed projective measurements, uncharacterized drifts in the detectors' effi- 
ciency, or a non-identical ensemble of unknown quantum states. The "precision" of a tomo- 
graphic reconstruction measures the statistical error in p, and is strongly dependent on the total 
number of measurable states N in the identical ensemble (which is in turn dependent on the en- 
tanglement source's pair production rate, R, and the total single-qubit measurement efficiency, 
7]). In general, the tomographic precision decreases as T (and therefore N) decreases fT3ll . For 
sufficiently small T we can neglect systematic effects and quantify tomographic precision (as a 




Fig. 1. (a) Tomographic precision F p (A^,p ideal ) for p ideal = |0 + )(0 + | with |0 + ) = 
-^(\HH) + | VV)). Each data point represents a 2000-tomography numerical Monte Carlo 
simulation of the average fidelity between the reconstructed density matrix and Pideab un_ 
der realistic assumptions about the system noise (a coincidence- to- accidental ratio of 3). 
Each simulated tomography utilizes four detectors and nine coincidence measurements, 
such that each unknown quantum state in the N- state ensemble is projected onto one of 
nine four-element orthonormal bases (e.g., HH, HV, VH, VV). Results for both the maxi- 
mum likelihood technique and the truncated-eignevalue, linear-least- squares-fit technique 
are shown. For a given N, the maximum likelihood technique is slightly more precise fT3l . 
(b) Using the same simulated data, F p is shown as a function of total tomography time T 
for two different experimental systems: a traditional free- space tomography system with 
77 = 0.1, z s = 5 s, % a = 5 s and an entangled photon polarimeter with 77 = 0.07, z s = 0.02 s, 
T a = 0.001 s. In both systems R = 10 6 pairs/second and M = 9. 



function of N and of Pideal) to be the average fidelity between p and Pideai 



F p (N, pideal) =F(p, Pideal) = ^Tr j yj y/p p idea i y/p jj . (2) 

Note that the equation above uses the usual definition for fidelity between two mixed states |[T4lL 
which for a pure Pid ea i = W){wV simplifies to the more familiar F (p,Pideai) = Tr{pp ide ai} = 
(y|p|V>- FigureUta) shows F p (N) for p ide ai = l4> + >(0 + | with |0+> = ^(\HH) + \W)) 9 where 
each data point represents a Monte Carlo simulation of the average fidelity between a recon- 
structed density matrix and the ideal unknown state. 

Two-qubit polarimetry is an application of two-qubit polarization tomography which maxi- 
mizes precision for very short T (< Is), allowing an experimenter to manipulate an entangled 
photon source using real-time tomographic feedback (by updating after every measurement, the 
time between updates can be reduced to T/9). (In this paper, entangled photon polarimetery 
refers to the application of two-qubit polarimetry to entangled photon states.) Because maxi- 
mizing precision requires maximizing N, the ideal entangled photon polarimeter will minimize 
both the time between measurements (z s ) and the time for numerical analysis (r a ): 

N = Rri 2 MT m = Rri 2 (T - Mx s - x a ) . (3) 

Although Eq. © can be used to derive the total time necessary to perform a single tomogra- 
phy with a given precision, an entangled photon polarimeter will likely perform many tomogra- 



phies in series. In this configuration, the tomographic measurements and the numerical analysis 
of those measurements can be parallelized in one of two ways. For x a < M(i m + z s ), a complete 
set of M measurements can be analyzed at the same time the next set of M measurements are 
being performed, leading to one tomographic result being displayed to the experimenter every 
M(t m + Ty) seconds. For T a < t m + z s , a tomographic result can be analyzed and displayed after 
every measurement, rather than after every complete set of M measurements. In other words, 
after every measurement, the previous M measurements are used to reconstruct an updated den- 
sity matrix, leading to a faster refresh rate based on a tomographic "rolling average". Similarly, 
this configuration can be altered in real time to utilize even more measurements (e.g., AM) for 
increased precision (analagous to averaging multiple traces on an oscilloscope). 

3. Experimental Details 

The entangled photon polarimeter presented here is based on a previous apparatus for free-space 
telecom-band quantum state tomography Il0l[T2|, which although accurate, is too slow to pro- 
vide real-time feedback. Three key improvements have dramatically improved the tomographic 
speed while maintaining precision: bulk wave plates have been replaced with fast electro-optic 
modulators (EOMs), an array of four single-photon detectors triggered at 8 MHz have been 
replaced with an array that is triggered at 50 MHz, and the traditional maximum likelihood 
reconstruction technique has been replaced with a much-faster linear-least-squares-fit method. 

Below, we briefly discuss the differences between these two techniques after reviewing the 
entangled photon source used to test the tomography apparatuses. Figure [Hb) highlights the 
differences between the two techniques, showing the expected tomographic precision F p as a 
function of total tomography time T. 

3.1. Entangled Photon Source 

To test the entangled photon polarimeter, we utilize a fiber-based, frequency-degenerate, 1550- 
nm, polarization entangled photon-pair source |[T0l . The source utilizes spontaneous four- wave- 
mixing in dispersion- shifted fiber and is pumped by 50-MHz repetition rate dual-frequency 
pulses spectrally carved from the output of a femtosecond pulsed laser. Because the output 
photons are identical, reverse Hong-Ou-Mandel interference in a Sagnac loop is used to deter- 
ministically split the output photons into separate output single-mode fibers. See Fig. 0a). 

The same source is used to test two separate tomography systems, the automated wave-plate- 
based apparatus first described in |[T0l and the entangled photon polarimeter presented here. 



Circ: Circulator EDFA: Erbium-doped fiber amplifier 

DSF: Dispersion shifted fiber FPC: Fiber polarization controller 

DGF: Double grating filter PBS: Polarizing beam splitter 

HWP: Half-wave plate PDD: Polarization dependent delay 




Fig. 2. (a) The entangled photon source used to test the entangled photon polarimeter. 
(b) The entangled photon polarimeter, composed of fast electro-optic modulators, in-fiber 
polarizers, and a four-detector array. 



3.2. Polarization Measurements 

Traditionally, two-qubit polarization tomography is performed using bulk, free-space, birefrin- 
gent crystals (i.e., wave plates). A quarter- and a half- wave plate followed by a polarizer on each 
of the two qubits can implement an arbitrary projective measurement lH2l . By collecting pho- 
tons from both the transmitted and the reflected ports of each qubit's polarizer, one can project 
an unknown photon pair into one of four orthonormal basis states, defined by the wave plates. 
If well characterized, this can lead to a very accurate tomography, though the measurement-to- 
measurement transition time T s will in general be very large (^ 5s). For the fiber-based source 
above, this type of polarization analyzer will lead to a single-qubit loss of « 1.5 dB (including 
the fiber to free-space to fiber coupling losses). 

To decrease z s , we have constructed an all-fiber/waveguide polarization analyzer based on 
electro-optic modulators (EOMs). These LiNb0 3 EOMs (EOSpace, model PC-B4-00-SFU- 
SFU-UL) allow precise control of both the retardance and optic axis of a birefringent crys- 
talline waveguide using the fringe fields from three electrodes. In general, this process has an 
extremely short response time leading to EOM switching rates of up to 10 MHz. In practice, 
we are able to implement arbitrary polarization measurements at 125 kHz, which is a limit set 
by the speed of our computer-controlled voltage sources. 

Although high-speed, EOMs are more difficult to precisely characterize than bulk wave 
plates; using a standard polarimeter we have characterized the six transformations performed 
by each EOM-based analyzer (corresponding to projections onto the H, V, D = (H + V)/ \/2, 
A = (H - V) /y/2, R = (H + /V) /y/2, and L = (H - /V) /y/2 basis states). EOM projections 
deviated from an ideal measurement by an average of 2.1 degrees on the Poincare sphere. The 
single-qubit losses of the EOM-based analyzers varied between 3.0-3.4 dB. 

3.3. Single-Photon Detection 

Single-photon detection is performed using a four-detector array of InGaAs avalanche photo- 
diodes (APDs) operated in the gated Geiger mode. By increasing the speed of these detectors 
from 8.3 MHz to 50 MHz, the entangled photon polarimeter achieves a 6-fold speed increase 
relative to previous implementations of quantum state tomography which utilized the same 
telecom-band detection systems. Moreover, by synchronizing the detector-array's control soft- 
ware with the EOM-based analyzers, we have reduced the switching time to z s = 20 ms. By 
upgrading the detector control software to eliminate extraneous electronic delays, we anticipate 
that this will approach the 125 kHz limit (t s = 10 /is) imposed by the EOM voltage controllers. 
The quantum efficiency of each detector at 1550-nm is approximately 20%, with a measured 
dark-count rate of 1-4 x 10 -4 per pulse. 

3.4. Tomographic Reconstruction 

Traditionally, the maximum likelihood technique has been used to reconstruct a two-photon 
state's density matrix from a series of coincidence measurements, which numerically solves 
for the density matrix p most likely to reproduce the measured counts (TU [121 . This method 
always produces a legal state, but is relatively slow (x a « 5 s). 

By using a simpler analysis technique based on a linear least-squares fit, we are able to 
increase the state reconstruction speed by more than three orders of magnitude fT3l . We use the 
2-qubit Stokes vector as a linear model, and solve the following least-squares problem: 

wM-S = wC. (4) 

Here, M is the set of measurements, which can be arbitrary POVMs; C is the measured counts, 
and S is the Stokes vector we solve for; w is a weight vector representing the distribution 
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Fig. 3. Selected frames from the nine fps video of a two-qubit photon state's transition 
from separability to entanglement. Each frame shows a density matrix reconstructed using 
the previous 36 measurements 4 s of data). 

width for each measurement. We assume the counting process to be Poissonian, and use the 
large-N limit where the Poisson distribution is approximated as a Gaussian with width \/N. To 
guarantee a legal density matrix, we post-process the least-squares fit by truncating the negative 
eigenvalues lTT3l . We have found that this type of linear fit provides results identical to those 
obtained via the maximum likelihood method with a negligible drop in precision (see Fig. [Ha)), 
only much faster (~ 1.3 ms per tomography using Matlab on a 2.4-GHz CPU). 

This three-order-of-magnitude speed increase allows us to display a new frame (i.e., to- 
mography result) after every measurement, reconstructed using the previous M measurements. 
For four-detector, complete-basis polarization analyzers (described above), only nine measure- 
ments are needed to perform a complete tomography. Note that it is often experimentally opti- 
mal to perform a redundant set of 36 measurements in order to detect and/or correct for system- 
atic errors such as source intensity drift, detector efficiency drift, or polarizer crosstalk fT2l . 

4. Entangled Photon Polarimeter Performance 

By utilizing fast EOM-based analyzers, a four-detector array triggered at 50 MHz, and a linear 
least-squares algorithm for tomographic reconstruction, the entangled photon polarimeter is 
capable of performing nine tomographies per second. Operated at this speed, T m = 80 ms, 
T s = 20 ms, and z a = l ms. Total single-qubit insertion loss is measured to be 7] = 3-3.4 dB 
(not including detector inefficiency). The tomographic precision is estimated using a Monte 
Carlo simulation of this polarimeter's application to the entanglement source pictured in Fig. [2] 
(resulting in w 1000 coincidences per second). For nine-measurement tomographies (T « 1 s), 
/>(p,Pideai) ~ 92%. For 3 6-measurement tomographies (T w 4 s), /^(p,Pideai) ~ 96%. 

To experimentally verify this performance, we recorded three 9-fps live videos of a two- 
photon polarization state using the 3 6-measurement configuration. First we recorded two videos 
where the measured state is not changed during the course of the measurement run, for a totally 
separable pure state, \HH), and a maximally entangled state, |0 + ). By analyzing each frame and 
comparing it to the target state, we directly measured the system precision to be 98% ±1% (for 
\HH}) and 95% ± 2% (for |0 + )). Finally, we recorded a video of a two-photon state's transition 
from separability to entanglement (the transition is physically implemented by rotating the wave 
plate HWP in the entangled photon source setup — see Fig.0. Selected frames from this video 
are shown in Fig. [3] 
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